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ABSTRACT 

We present velocity-delay maps for optical H I, He I, and He II recombination lines in Arp 151, recovered 
by fitting a reverberation model to spectrophotometric monitoring data using the maximum-entropy method. 
H I response is detected over the range 0-15 days, with the response confined within the virial envelope. The 
B aimer-line maps have similar morphologies but exhibit radial stratification, with progressively longer delays 
for H7 to H/3 to Ha. The He I and He II response is confined within 1-2 days. There is a deficit of prompt 
response in the B aimer- line cores but strong prompt response in the red wings. Comparison with simple models 
identifies two classes that reproduce these features: freefalling gas, and a half-illuminated disk with a hotspot 
at small radius on the receding lune. Symmetrically illuminated models with gas orbiting in an inclined disk 
or an isotropic distribution of randomly inclined circular orbits can reproduce the virial structure but not the 
observed asymmetry. Radial outflows are also largely ruled out by the observed asymmetry. A warped-disk 
geometry provides a physically plausible mechanism for the asymmetric illumination and hotspot features. 
Simple estimates show that a disk in the broad-line region of Arp 151 could be unstable to warping induced by 
radiation pressure. Our results demonstrate the potential power of detailed modeling combined with monitoring 
campaigns at higher cadence to characterize the gas kinematics and physical processes that give rise to the broad 
emission lines in active galactic nuclei. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individual (Arp 151) 



1. INTRODUCTION 

Active galactic nuclei (AGNs) have long been known 
to undergo dramatic and rapid variations in bri ghtness 
(iMatthews & Sandagd [T963I: iSmith & H"offleE [T963h . This 
variability is exploited with the technique o f reve rberation 
mapping (Blandford & McKee 1982; Peterson 1993h to study 
the broad-line region (BLR) of AGNs, a region of gas only 
~ 0.01 pc in size and presumed to be photoionized by the 
emission from an accretion disk around a supermassive black 
hole. Spatially unresolvable even for the most nearby galax- 
ies, the BLR is studied in the time domain through monitor- 
ing the variable continuum emission (presumably originating 
close to the black hole) and its reprocessing in the BLR gas, 
discernible as variations in the broad emission-line fluxes and 
profiles. 

Under the assumptions that (1) the continuum-emission re- 
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gion is much smaller than the BLR, (2) time delays arise pre- 
dominantly from light- travel time, and (3) the relationship 
between the ionizing and observed continuum is smooth, a 
linearized echo model can be used to express changes in the 
emission lines (AL) as time-delayed responses to changes in 
the continuum (AC): 



AL(v,0 
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Here, AL(v, t) is the change in the emission-line flux at line- 
of- sight velocity v and time t, AC(t) is the change in the con- 
tinuum flux at time t, and ^(v, r) is the so-called transfer func- 
tion, or velocity-delay map, giving the distribution of the line 
response over line-of- sight velocity v and time delay r (i.e., 
the emission-line response to a (5-function continuum event). 
Information on the BLR geometry, kinematics, and ionization 
structure is encoded in ^(v,r), making recovery of velocity- 
delay maps a key goal of echo-mapping experiments. 

Due to the formidable requirements of long-duration 
monitoring, high temporal sampling, and high-quality and 
homogeneous data for recovery of velocity-delay maps 
(JHorne et al.ll2QQ4|) . most echo-mapping experiments have re- 
stricted the analysis to a simpler problem: measuring the 
mean delay r, and hence the mean radius R ^ cr of the 
BLR, through cross-correlation of the line and continuum 
light curves (e.g.. lAntonucci & Cohen 1981; iPeterson et al.l 
11983"; "Kaspi et alJ|2000;|Peterson et al. 2004). Analyses com- 
paring delays in different parts of the emissio n-line profiles , 
such as the red wing versus t he blue wing (e.g.,lGaskelllll988l: 
iKoratkar & GaskelL 119891: ICrenshaw & Blackwelll 119901) . 
have sometimes shown evidence of asymmetry, hinting at 
gross characteristics (inflow vs. outflow vs. circulation) of 
the BLR gas flow. 
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For certain high-quality data sets, delay maps (^(r) = 
/^(v,r)(iv) have been recovered. These include delay maps 
for hya and C IV, among other lines, derived from the In- 
ternational Ultraviolet Explo rer (lUE) monitoring dataset for 
NGC5548 (Krolik et al. 1991), and the delay maps recovered 
by lHorne et al. (1991) from ground-based optical monitoring 
of H/3 in NGC 5548. The delay maps generally agree with the 
results of cross-correlation analysis, while allowing add itional 
insights into the details of the BLR. lHorne et al.l (Il991h find a 
lack of H/3 response at zero time delay in NGC 5548, a signa- 
ture of anisotropic (inward) H/3 response and/or a paucity of 
BLR gas near the line of sight. 

In addition to prese nting delay maps for se veral ultravio- 
let lines in NGC 4151. luirich & Horng (1 1 9961) also present a 
partially recovered velocity-delay map for C IV, which un- 
fortunately suffered from strong intrinsic absorption in the 
line core throughout the campaign. Three separate groups 
attempted to recover velocity-delay maps from an lUE and 
Hubble Space Telescope monitoring campaign of NGC 5548; 
however, these studies were again hampered by the quality 
of the C IV data us ed, and the final c onclusions ranged from 
no radial motion (Wanders et al."1995) to some radial infall 
(^one & Krolik 1996) to radial outflow (Chiang & Murrav 
[19961: iBottorff et al.l[T997h . 

Recent results from the Lick AGN Monitoring Project 
(LAMP; Bentz et al. 2008, 2009) and from MDM Obser- 
vatory (Penne y et al.i i2009a.b) have shown that recovery of 
velocity-delay maps may be within the reach of current AGN 
monitoring programs due to their monitoring baselines of 
several months, daily sampling, and homogeneous data with 
high signal-to-noise ratios (S/N). In this Letter, we present 
velocity-delay maps recovered with maximum-entropy tech- 
niques for six optical H and He recombination lines in the 
LAMP spectra of Arpl51 (Mrk40): Ha, H/3, H7, He I 
A5876, He I A4471, and He II A4686. We also suggest plau- 
sible models for the geometry and kinematics of the BLR in 
Arpl51 that could account for the observed velocity-delay 
structure. 

2. OBSERVATIONS 

Details of the p hotometric light curves are presented by 
IWalsh et al.l (120091) . In brief, broad-band Johnson B and V 
images of Arp 151 were obtained at the 0.8-m Tenagra II tele- 
scope in southern Arizona between 2008 February 27 and 
May 16 (UT here and throughout). Standard reduction tech- 
niques were employed, and photometric light curves were 
measured with differential photometry relative to fi eld stars. 
Absolute flux calibrations were determined using iLandoltl 
|l992) standard- star fields. 

Spectroscopic monitoring of Arp 151 was carried out at 
the Lick Observatory 3-m Shane telescope with the Kast 
dual spectrograph between 2008 March 25 and May 21 (see 
lBentzetal.i i2009 for details). The Kast red-side CCD with 
the 600 lines mm~^ grating allowed spectral coverage over 
the range 4300-7100 A. Spectra were obtained at a fixed po- 
sition angle of 90° through a 4''-wide slit, with a typical S/N 
^100 per pixel at rest- frame 5100 A. Flux calibrations were 
determined from nightly spectra of standard stars. 

Small time- variable and wavelength-dependent corrections 
were made to the spectra to account for three types of sys- 
tematic error: (1) wavelength shifts due to flexure and atmo- 
spheric differential refraction, (2) spectral blurring due to at- 
mospheric seeing and instrumental resolution, and (3) photo- 
metric errors due to slit losses and atmospheric transmission 



(see Bentz et al. 2010, in preparation, for a detailed descrip- 
tion of the methodology). After correcting for these errors, the 
continuum- subtracted line-profile variations were subjected 
to echo-mapping analysis with MEMECHO. 

3. MEMECHO LIGHT-CURVE MODELING 



MEMECHO rsee lHorneet"an[T99Tl: iHornd [T9941 for de- 
tails) uses a maximum-entropy technique to recover ^(v,r) 
by fitting the linearized echo model of Equation ([T]) to the ob- 
served continuum light curve and responding emission-line 
velocity profiles. MEMECHO finds the "simplest" positive 
image pt that fits the data. The image pi includes an evenly 
sampled driving continuum light curve C(t), a background 
spectrum Lo(A) representing nonvariable line emission, and 
the velocity-delay map ^(v,r). The quality of the fit to the N 
data points is measured by x^ summed over measurements 
of the driving continuum light curve and the reverberating 
emission-line spectrum. Simplicity is measured by the en- 
tropy S = ^iPi-qi-piln(pi/qi) relative to a "default image" 
qi. Since dS/dpt = -^-ln(pi/qi), maximum entropy at 5" = 
occurs when pt = qt. Our default image qt slightly blurs the 
image pt (for example qt = .^/pJITPi+i for a one-dimensional 
image), and thus S increases toward for smooth positive 
images. With suitable datasets MEMECHO delivers a well- 
defined map maximizing S for a specified x^ /N. The reso- 
lution of the map and the quality of the fit improve as x^ /N 
is reduced. When x^ is set too low, the fit attempts to follow 
noise in the data, resulting in noisy maps. The maps we show 
are for x^/N = 1.1, giving the highest resolution maps and 
the best fits we could achieve while avoiding the introduction 
of spurious structure, typically large excursions in the gaps 
between measurements of C(t). 

We chose the 5-band photometric light curve as the driving 
continuum light curve for the MEMECHO analysis. The vari- 
ations in B are stronger t han in V with no evident time delay 
between the two (IWalsh et al. 2009), and the photometric light 
curves are more accurately calibrated than the spectroscopic 
continuum light curve, which suffers from aperture effects. 
Furthermore, the contribution of the broad e mission line s to 
the photometric variability is negligible (see IWalsh et aP for 
a complete discussion). The V-band light curve and the in- 
tegrated Ha, H/3, and Hj light curves were then modeled as 
echoes of B, as were the emission-line light curves at each 
wavelength. Thus, we fit delay distributions simultaneously 
to a total of 1373 individual light curves — one per spectral 
pixel (each with AA = 2 A), plus the V-band and integrated 
B aimer-line light curves. We enforce causality by requiring 
non-negative delays, Tmin = Odays. We set Tmax = 20 days, 
which is ~ 1/3 of the spectroscopic monitoring baseline, and 
we use a uniform spacing of Ar = 0.5 day in the delay maps 
and continuum light curve. The linearized echo model gives a 
tangent-line approximation to what may in principle be non- 
linear responses to changes AC(t) = C(t) - Cq away from a 
reference continuum level Q. For Q we adopt the mean of 
the observed 5-band fluxes. 

In Figure 1, we show selected light curves of Arp 151 along 
with their MEMECHO fits and recovered delay maps. The 
selected wavelengths sample the cores of the emission lines 
plus the B aimer-line wings. The 5-band (driving) light curve 
is shown in the bottom panel. The value of x^/^ is reported 
for each light curve shown, and is 1 . 1 when summed over all 
light curves including those not shown. 
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4. THE VELOCITY-DELAY MAPS 

The delay maps at each wavelength are assembled into a 
velocity-delay map shown over the full spectral range in Fig- 
ure 2 and centered on the six strongest emission lines in Fig- 
ure 3. To aid in comparing the velocity-delay structure of dif- 
ferent emission lines, Figure 4 presents false-color maps with 
red, green, and blue colors giving the velocity-delay distribu- 
tions of different emission lines as indicated. 

Taken together, the velocity-delay maps exhibit a virial 
structure with the emission-line response covering a wider 
range of velocities at smaller delays. This is perhaps best 
seen in Figures 3 and 4, where the line responses are con- 
fined within the "virial envelope" v^ = GM/cr for Myinai = 
1.2 X lO^M© (Bentzetal. 2009). 

The B aimer-line maps all have the same basic structure, 
with longer delays in the core and smaller delays in the line 
wings, and with the red-wing response much stronger than 
that of the blue wing. The Ha response is detected from 
days in the red wing to about 15 days in the line core. The 
other B aimer lines have a similar structure, but with H/3 ex- 
tending to about 7 days and H7 to 5 days. The mean de- 
lays from these maps are consistent with the mean lag times 
from the cross-correlation analysis presented by B entz et al.l 
(2010). A plausible interpretation of this effect is a radial 
stratification resulting from optical-depth effects within the 
B aimer series ( Kori sta & Go ad 2004). 

The response of the helium lines is consistent with zero 
time delay, except for some response at small positive de- 
lay sjnJhemapofHe I A5876 and He I A4471. As discussed 
by 'Bentz et al.' (^2010), the helium response occurs on shorter 
timescales than our monitoring sampling and is therefore un- 
resolved. 

The Ha map shows a "curl" at long time delays on the blue 
side of the line outside the virial envelope. This weak feature 
may be an artifact, as it does not always appear in maps made 
with slightly different fitting parameters. Features that appear 
in all of the Balmer lines, such as the strong emission in the 
red wings, are more secure. 

The maps show a lack of prompt response in the core of 
all three Balmer lines. The delay map for the integrated H/3 
profile in Figure 1 shows a paucity of re sponse at zero lag , 
similar to what was found for NGC 5548 (iHorneet al.ll 19911) . 
even though the full velocity-delay map clearly shows gas re- 
sponding with zero lag in the red wing of H/3. 

5. DISCUSSION 

We have examined the predicted velocity-delay structure 
for a variety of simple BLR models to identify classes of 
models that fail and succeed in qualitatively reproducing the 
main features of the observed velocity-delay maps. Two suc- 
cessful classes of models are represented by the Freefall and 
Disk + Hotspot models shown in Figure 5. For all four models 
in Figure 5 we ad opt Mbh = 7 x 10^ M© (assuming / = 5.5; 
lBentzetal.l[2009h . 

A thin spherical shell of radius R infalling with veloc- 
ity V covers a sloped line on the velocity-delay map, with 
r = {l-co^O)R/c and v = Vco^O, the delay increasing from 
r = at V = +y on the near side to r = 2R/c at v = -V on 
the far side of the shell. The Freefall model in Figure 5 has 
a spherical distribution of infalling gas with Ri^ = 0.1 It-days 
and 7?out = 4 It-days, and V = y^GM/R, giving small delays 
on the red side and longer delays with a virial envelope on 
the blue side of the map. This naturally produces the required 



asymmetry. The lack of prompt response in the line core can 
arise by means of an inward radiation anisotropy (as displayed 
in Figure 5 with Fin/^totai = 0.8) and/or by reducing the re- 
sponse of the inner gas. Models with radially outflowing gas 
have the opposite asymmetry and are thus largely ruled out. 

For an inclined circular orbit the velocity-delay structure is 
an ellipse with r = (1 + sin /cos 0)R/c and v = V sin / sin 0. The 
basic virial envelope, with a wider velocity range at smaller 
delays, can be reproduced by distributions of circular Ke- 
plerian orbits ranging from isotropically distributed angular 
momentum to inclined disks (e.g., the Thick Spherical Shell 
and Isotropically Illuminated Disk models in Figure 5). For 
such models the Balmer response at small delay in the line 
core can be reduced by (1) anisotropic ionizing radiation di- 
rected toward the far side, (2) anisotropic inwardly directed 
line response, and/or (3) excluding edge-on orbits to form an 
inclined disk-like geometry. The enhanced prompt response 
in the red wing is hard to reproduce with a symmetric distri- 
bution of circular orbits. If the orbits circulate in the same 
direction, however, then the required asymmetry can be in- 
troduced by azimuthal structure-enhancing response (a "hot 
spot") on the receding side of the inner disk. While this is 
somewhat ad hoc, it might plausibly arise in the context of 
a warped-disk geometry, the warp exposing gas to ionizing 
radiation at small radius while shielding it at larger radius. 
The Disk + Hotspot model displayed in Figure 5 is inclined 
(/ = 20°) with 7?in = 0.1 It-days and 7?out = 4 It-days. No emis- 
sion originates from the near side of the disk due to its warped 
shape, and excess emission arises on the receding portion of 
the inner disk at 0.1-0.3 It-days. 

Both the Infall model and Disk + Hotspot model qualita- 
tively reproduce the observed features of the velocity-delay 
maps, and more work is necessary to determine which family 
of models is preferred. In the meantime, we note that warped 
disks are seen in many astrophysical environments and have 
previously been invoked to explain various AGN phenom- 
ena (e.g., flux va riations i n the double-peaked broad emission 
from NGC 1097. IStorchi- Bergmann et al. 1 9971: Fe emiss ion- 
line profiles in the X-rays, Hartnoll & Blackman"2000'; the 
misalignment of radio jets and nuclear disks, Schmitt e t al.l 
120021; and a sel f-contained apparatus for AGN unification, 
e.g., iNay akshini i200 5 ) . We c onsider here the plausibility of 
irradiation-induced warping ( Pringldll996h . in which an op- 
tically thick, geometrically thin disk irradiated by a central 
source is unstable to warping from radiation pressure. The 
instability criterion is 



R 



> ^TT^rj^e-^ 



(2) 



where Rq = GMbu/c^, t] is the ratio of the viscosities in the 
plane of and perpendicular to the disk (assumed here to be 
equal, i.e., r] = l), and e is the accretion efficiency. Assuming 
e ^ 0. 1 , a disk in the BLR of Arp 151 could become unstable 
on length scales R/Rg > 800 ^ 0.3 It-day, inside the scale 
of the optical BLR. For typical estimates o f BLR parame- 
ters (following IStorchi-Bergmann et al.|[T997h . the precession 
timescale would likely be decades to centuries. 

6. SUMMARY 

We have recovered velocity-delay maps for six optical hy- 
drogen and helium recombination lines in the LAMP spec- 
tra of Arp 151 using the maximum-entropy method. These 
are the most detailed velocity-delay maps constructed to date. 
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The individual B aimer lines have similar velocity-delay struc- 
ture, and the maps show responses with delays increasing 
from H7 to H/3 to Ha. This can be interpreted as radial 
stratification through the B aimer series and is expected from 
optical-depth effects. All three lines show prompt response in 
their red wings, contrary to what has been seen in NGC 5548. 
The helium-line response is mostly unresolved in time. 

The features seen in the B aimer-line velocity-delay maps 
can be reproduced qualitatively with either of two simple ad 
hoc models: a freef ailing BLR, or a partially illuminated thin 
disk with a localized excess of emission at small radii that 
could be interpreted as a warped disk. Warped disks are 
known to occur in many environments and on many astro- 
nomical scales, and it is plausible that the BLR in AGNs is 
another manifestation of this common configuration, as has 
been previously suggested in the literature. It appears to be 
impossible for a simple disk model or a thick spherical shell 
composed of randomly inclined circular orbits to accurately 
reproduce the features seen in the velocity-delay maps, in- 
cluding the strong emission at short lag times in the red wings 
of the B aimer lines and the lack of emission at short lag times 
in the line cores. Radial outflows are also largely ruled out by 
the observed asymmetric response. 

The simple models and plausibility arguments presented 
here are the first steps in interpreting the velocity-delay maps 
for Arp 15 L We are currently investigating other LAMP tar- 



gets and expect to recover velocity-delay maps for a few ad- 
ditional objects, which may lead to insights into BLR dif- 
ferences across the Seyfert 1 population. Future monitoring 
programs will benefit from even higher temporal sampling 
and the avoidance of any gaps in the driving continuum light 
curve, allowing the driving light curve to provide stronger 
constraints and leading to more detailed velocity-delay maps 
of the emission-line responses. More sophisticated model- 
ing will also help with interpreting the velocity-delay maps 
and determining which families of BLR models can be ruled 
out, allowing additional insights into the complexities of AGN 
BLRs. 
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Fig. 1 . — MEMECHO fits to the light curves at selected wavelengths in the spectra of Arp 151 {right panel) and their delay maps (left panel). The rest-frame 
central wavelength of each pixel is listed with its delay map, and the x^ /N values are listed with the light curves. We include four supplementary integrated light 
curves (V band, Ho;, H^, and H7) as fitting constraints, which are displayed near the bottom. Multiple measurements in the same night are averaged together in 
the integrated light curves, but are treated separately for the non-integrated light curves. 
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Fig. 2. — Velocity-delay map ^(A,r) over the full wavelength range (greyscale), and projections onto the wavelength axis ^(A) = f'^(X,T)dT (bottom) and 
the time-delay axis ^(r) = / ^(A, T)dX (right). 
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Fig. 3. — Velocity-delay maps for each of the six optical H and He recombination lines in the LAMP spectra of Arp 151. The dotted lines in each panel show 
the "virial envelope" V^rc/G = 1.2 x IO^Mq, based on the "virial product" of time lag and line width for H/3 jBentz et d1l2009l) . 
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Fig. 4. — Left: Comparison of the velocity-delay maps for Ha (red), H^ (green), and H7 (blue). Right: Comparison of the velocity-delay maps for H/3 (red), 
He I A5876 (green), and He II A4686 (blue). As in Figure 3, the dotted lines in each panel show the virial envelope. The Balmer lines all have the same basic 
shape, but are graduated in their response ranges. The helium-line response is mostly unresolved in time. 
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Fig. 5. — Velocity-delay map for H^ compared to simple models for Mbh = 7 x 1O^M0: (top) a freefalling BLR with inward emission anisotropy (left), and 
a partially illuminated disk with a hot spot at small radii on the receding side (right); (bottom) a biconically illuminated thick shell with circular orbits in a 
Keplerian potential and inward emission anisotropy (left), and a fully illuminated disk with inward emission anisotropy (right). The top models more accurately 
reproduce the features seen in the H/3 velocity-delay map than do the bottom models. All models have been smoothed to match the lower resolution of the 
recovered velocity-delay maps. The dotted lines in the H^ panel again show the virial envelope, and the red data points show the average lag measurement per 
velocity bin from Bentz et al. ( 2009) . 



